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We use the recently proposed jettiness-subtraction scheme to provide the complete calculation 
of Higgs boson production in association with a jet in hadronic collisions through next-to-next- 
to-leading order in perturbative QCD. This method exploits the observation that the A-jettiness 
event-shape variable completely describes the singularity structure of QCD when final-state colored 
particles are present. Our results are in agreement with a recent computation of the gg and qg 
partonic initial states based on sector-improved residue subtraction. We present phenomenological 
results for both fiducial cross sections and distributions at the LHC. 


INTRODUCTION 

The continued investigation of the Higgs boson will be 
a top priority of the high energy physics community dur¬ 
ing the coming years. The comparison of Standard Model 
(SM) predictions with data from Run I of the Large 
Hadron Collider (LHC) was limited by the statistical pre¬ 
cision of the experimental data. This will no longer be 
the case during Run II, and systematic errors will dom¬ 
inate. The largest systematic error currently hindering 
our understanding of Higgs properties is the theoretical 
understanding of the SM prediction. This is the case for 
the well-measured di-boson modes [1] which dominate 
the overall signal-strength determination. The theoreti¬ 
cal uncertainties must be reduced in order to sharpen our 
understanding of the mechanism of electroweak symme¬ 
try breaking in Nature. 

Improvements in both the overall production rate of 
the Higgs boson and in the modeling of its kinematic 
distributions are needed to match the expected experi¬ 
mental precision of Run II. The total cross section is cur¬ 
rently known at next-to-next-to-leading order (NNLO) 
in the strong coupling constant [2-4], and the resumma¬ 
tion of enhanced threshold logarithms at higher orders 
is known [5-7]. The N 3 LO corrections have recently be¬ 
come available [8], as are phenomenological predictions 
incorporating much of our current knowledge of the inclu¬ 
sive cross section [9] . To improve the modeling of Higgs 
kinematics at the LHC, precision predictions for Higgs 
production in association with jets are needed. In decay 
modes such as ff —> W + W~ —> l + vl~v in which a mass 
peak cannot be reconstructed, the need for theory is crit¬ 
ical. The resummation of jet-binning logarithms beyond 
NLO in the WW final state has been shown to lead to a 
factor of two reduction in the theory uncertainty affecting 
this channel [10], indicating that higher-order corrections 
play a vital role in the analysis of Higgs properties in this 
channel. 

Initial results for Higgs production in association with 


a jet at NNLO are available, including both partial re¬ 
sults in pure Yang-Mills [11, 12] and a recent compu¬ 
tation of the complete gg and qg partonic scattering 
channels [13]. We present here a full computation of all 
partonic channels through NNLO for the Higgs plus jet 
process using the recently proposed jettiness subtraction 
scheme [14]. Our aim in this manuscript is two-fold; first, 
we wish to demonstrate the application of jettiness sub¬ 
traction in a non-trivial example. Higgs production in 
association with a jet provides such a test case. Second, 
given the importance of the precision Higgs program at 
the LHC, we wish to confirm the recent computations of 
this process. We find agreement with the results obtained 
by the calculation in Ref. [13] . Our result in addition con¬ 
tains the quark-initiated scattering processes at NNLO, 
which were not included in previous calculations [13]. We 
confirm that they have minimal phenomenological im¬ 
pact. 

Our paper is organized as follows. We begin by review¬ 
ing the recently introduced jettiness-subtraction formal¬ 
ism. The checks of our calculation, both internal ones 
and those involving other known results in the literature, 
are discussed next. We then present numerical results 
for fiducial cross sections and for several distributions of 
both the Higgs boson and the leading jet. We conclude 
with a summary and a discussion of potential future ex¬ 
tensions of the work presented here. 


REVIEW OF JETTINESS SUBTRACTION 

We review here the salient features of the jettiness- 
subtraction scheme for NNLO calculations, which was 
recently introduced in the context of the NNLO com¬ 
putation of W-boson production in association with a 
jet [14]. We begin by defining the A-jettiness event shape 
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variable Tn-, first introduced in Ref. [15]: 



(1) 


The subscript N denotes the number of jets in the final 
state. For the H +jet process considered here, we have 
N = 1. Values of 7i near zero indicate a final state con¬ 
taining a single narrow energy deposition, while larger 
values denote a final state containing two or more well- 
separated energy depositions. The first of these configu¬ 
rations will eventually reconstruct to one jet after impos¬ 
ing a jet algorithm, while the second will reconstruct to 
two or more jets. The pt are light-like vectors for each of 
the initial beams and final-state jets in the problem, while 
the qk denote the four-momenta of any final-state radia¬ 
tion. The Qi characterize the hardness of the beam-jets 
and final-state jets. We set Qi = 2 Ei, twice the energy 
of each jet. 

The cross section for Tn less than some value T§ ut can 
be expressed in the form [16, 17] 


cr(T N <TN Ut ) = J H®B®B®S® 


N 


IP- 
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The function H contains the virtual corrections to the 
process. The beam function B encodes the effect of radi¬ 
ation collinear to one of the two initial beams. It can be 
written as a perturbative matching coefficient convoluted 
with a parton distribution function. S describes the soft 
radiation, while J n contains the radiation collinear to 
a final-state jet. The ellipsis denotes power-suppressed 
terms which become negligible for Tn Qi- Each 
of these functions obeys a renormalization-group equa¬ 
tion that allows logarithms of Tn to be resummed. If 
this expression is instead expanded to fixed-order in the 
strong coupling constant, it reproduces the cross section 
for low Tn- The derivation of this factorization theo¬ 
rem in the small-7Ar limit relies upon the machinery of 
Soft-Collinear Effective Theory [18]. 

The basic idea behind jettiness subtraction is that Tn 
fully captures the singularity structure of QCD ampli¬ 
tudes with final-state partons. This allows us to calcu¬ 
late the NNLO corrections to processes such as Higgs plus 
jet in the following way. We divide the phase space ac¬ 
cording to whether Tn is greater than or less than Tff 1 * ■ 
For Tn > T^ ut there are at least two hard partons in 
the final state, since all singularities are controlled by 
jettiness. This region of phase space can therefore be 
obtained from a NLO calculation of Higgs production in 
association with two jets. Below T§ ut , the cross section 
is given by the factorization theorem of Eq. (2) expanded 
to second order in the strong coupling constant. As long 
as T^ ut is smaller than any other kinematic invariant 
in the problem, power corrections below the cutoff are 
unimportant. 


All ingredients of Eq. (2) are known to the appropriate 
order to describe the low Tn region through second or¬ 
der in the strong coupling constant. The two-loop virtual 
corrections are known for Higgs plus jet [19]. The beam 
functions are known through NNLO [20, 21], as are the 
jet functions [22, 23] and soft function [24], It is there¬ 
fore possible to combine this information to provide the 
full NNLO calculation of Higgs production in association 
with a jet. 


VALIDATION OF THE CALCULATION 

We describe here the various checks we have performed 
on our calculation. For the phase space region above 
Tx ut we need an NLO calculation of Tf+2-jets. We use 
MCFM [25, 26] for this purpose. Below T^ ut , we re¬ 
quire several separate terms. We have checked our im¬ 
plementation of the two-loop virtual corrections against 
those contained in PeTeR [27]. Our comparison of the 
soft function against known results in the literature is 
detailed in Ref. [24], 

An important check of our formalism is the indepen¬ 
dence of our final result from T{' ut ■ after both phase-space 
regions have been summed. This requires a consistent 
implementation of the jet, beam, and soft functions be¬ 
low the cut, and the NLO result for 77+2-jets above the 
cut. This also allows us to define the appropriate range of 
T cut for which power corrections are negligible. We show 
in Fig. 1 the NNLO correction to the cross section as a 
function of Ti ut ■ We also show separately the contribu¬ 
tions above and below Tf ut - The cross sections from the 
separate regions are each a factor of several larger than 
their sum, and vary by more than a factor of two over 
the Ti ut range studied. The sum is extremely stable over 
the entire range studied. 


NUMERICAL RESULTS 


We now present numerical results for Higgs production 
in association with a jet. We focus on 8 TeV proton- 
proton collisions in this paper. Jets are reconstructed 
using the anti-fc-r algorithm [28] with R = 0.5. We 
show results using the NNPDF [29] parton distribution 
functions. We use the perturbative order of the PDFs 
that is consistent with the partonic cross section under 
consideration: LO PDFs with LO partonic cross sec¬ 
tions, NLO PDFs with NLO partonic cross sections, and 
NNLO PDFs with NNLO partonic cross sections. We set 
the renormalization and factorization scales equal to the 
mass of the Higgs boson, pn = pp = mn- To estimate 
the residual theoretical error, we vary these scales simul¬ 
taneously around this central value by a factor of two. 
We set the mass of the Higgs boson as rriH = 125 GeV. 
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Figure 1: The separate cross sections for the regions 71 > 
T-Z ut and 7i < T\ ut , together with their sum, as a function 
of T\ ut , normalized to the NLO cross section. The lower 
panel shows the relative correction also with respect to the 
NLO cross section. The solid lines denote the result for the 
central scale choice p = ran, while the bands show the result 
as the scale as varied in the range tuh/2 < p < 2 m,H- The 
black dashed lines denote <j/(Jnlo = 1 (upper panel) and 
A a/aNLO = 0 (lower panel). 


We begin by studying the fiducial cross section for 
Higgs+jet production, which we define by imposing the 
following cuts on the final state jet: p^ > 30 GeV, 
\Y^ et \ < 2.5. Our results are shown in Table I. We 
can compare these results against the calculation of 
Ref. [13], which is based on the technique of sector- 
improved residue subtraction [30, 31]. The result quoted 
in Ref. [13] does not include a cut on |Y’ jet |, and also 
does not include the quark-initiated partonic channels 
qq , qq, and qq. Incorporating the |Y e *| cut in the sector- 
improved subtraction calculation [32] and removing the 
quark-initiated channels from our result, we find agree¬ 
ment within numerical errors. We note that the quark- 
initiated partonic channels reduce the NNLO result by 
approximately 1.5% within the fiducial region studied 
here, indicating that they have small phenomenological 
impact. 


p T et > 30 GeV, |F iet | < 2.5 

Leading order: 

3-llJj pb 

Next-to-leading order: 

4.8t 0 .9 Pb 

Next-to-next-to-leading order: 

5-5±g;2 pb 


Table I: Fiducial cross sections, defined by p^f* > 30 GeV, 
|W et | < 2.5, using NNPDF PDFs at each order of pertur¬ 
bation theory. The central scale choice is p = mrr- Results 
for p = mH/“i and p = 2rriH are shown as superscripts and 
subscripts, respectively. 

We next show several distributions in Higgs plus jet 


production. We show in Fig. 2 the rapidity distribution 
of the leading jet at each order in perturbation theory, 
as well as the A'-factors (defined as ratios of higher-order 
cross section over the lower-order ones) in the lower inset. 
The NLO corrections exhibit a slight shape dependence, 
with the corrections approximately 10-20% larger in the 
central region than near \Y^ et \ = 2.5. The NNLO correc¬ 
tions are flatter as a function of rapidity, and the NNLO 
distribution is entirely contained within the NLO scale 
variation band. In Fig. 3 we show the transverse momen¬ 
tum distribution of the leading jet. There is a shape de¬ 
pendence to the corrections, with the AT-factor decreasing 
as py f is increased. This trend is visible when going from 
LO to NLO in perturbation theory, and also when going 
from NLO to NNLO. We note that the NNLO result is 
entirely contained within the NLO scale-variation band. 
The shape dependence and magnitude of the NNLO cor¬ 
rections for the p^ et distribution are in agreement with 
the results of Ref. [13]. 



Figure 2: The rapidity of the leading jet at LO, NLO, and 
NNLO in the strong coupling constant. The lower inset shows 
the ratios of NLO over LO cross sections, and NNLO over 
NLO cross sections. Both shaded regions in the upper panel 
and the lower inset indicate the scale-variation errors. 

Finally, we show in Fig. 4 the transverse momentum 
of the Higgs boson in the ff+jet process. The NLO cor¬ 
rections range from 40% to 120% near pif = 60 GeV, 
depending on the scale choice. The magnitude of this 
correction decreases as the transverse momentum of the 
Higgs increases. The NNLO corrections are more mild, 
reaching only 20% at most for the central scale choice 
p = mu■ They also decrease slightly as the transverse 
momentum of the Higgs increases. The shape depen¬ 
dence and magnitude of the NNLO corrections for the p ^ 
distribution are in agreement with the results of Ref. [13]. 
We note that we have combined the two bins closest to 
the boundary pif = 30 GeV to avoid the well-known Su- 
dakov shoulder effect [33]. 
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Figure 3: The transverse momentum of the leading jet at LO, 
NLO, and NNLO in the strong coupling constant. The lower 
inset shows the ratios of NLO over LO cross sections, and 
NNLO over NLO cross sections. Both shaded regions in the 
upper panel and the lower inset indicate the scale-variation 
errors. 



Pt [ GeV] 

Figure 4: The transverse momentum of the Higgs boson at 
LO, NLO, and NNLO in the strong coupling constant. The 
lower inset shows the ratios of NLO over LO cross sections, 
and NNLO over NLO cross sections. Both shaded regions 
in the upper panel and the lower inset indicate the scale- 
variation errors. 


CONCLUSIONS 

We have presented in this manuscript a complete cal¬ 
culation of Higgs production in association with a jet 
through NNLO in perturbative QCD. Our computation 
uses the recently proposed method of jettiness subtrac¬ 
tion, a general technique for obtaining higher-order cor¬ 
rections to processes containing final-state jets. We con¬ 
firm and extend a recent calculation of the dominant 


gg and qg partonic channels through NNLO [11], and 
present additional phenomenological results for 8 TeV 
LHC collisions. We also present several distributions for 
the Higgs and the leading jet that can be measured with 
LHC data. Our results indicate that the perturbative se¬ 
ries is under good control after the inclusion of the NNLO 
corrections. We look forward to the comparison of our 
theoretical prediction with the upcoming data from Run 
II of the LHC. 
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